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Abstract 

Deep inelastic processes at small x are discussed in the framework of perturba- 
tive QCD at high energy. New results are presented on the quark anomalous 
dimensions beyond the leading logarithmic approximation, and their relevance 
to the structure functions being measured at HERA is pointed out. 
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1. Introduction 



Quantitative tests of QCD and searches for new physics at present and future hadron 
colliders are carried out at an increasingly large centre-of-mass energy y/S. In this regime 
a new kinematic region opens up, characterized by small values of the ratio x = Q 2 /S 
between the typical momentum \/Q^ transferred in the process and the total energy y/S. 
The most striking phenomenological feature at low x is the rise of the deep inelastic cross 
sections. This behaviour is qualitatively expected in QCD, and has been observed by the 
HERA collaborations Q, who have measured the structure function F 2 down to x- values 
as low as 1CT 3 -r- 1CT 4 . The observed rise is indeed steeper than any rise one finds in total 
cross sections for soft hadronic processes. Therefore it can only be accounted for by a hard 
QCD component, i.e., QCD interactions occurring at distance scales much smaller than 
Aq 1 cd . Such scales are controlled by perturbation theory. Let us then briefly recall the 
main features of QCD perturbation theory in the region of small x. 

Perturbative QCD predictions rely upon the factorization theorem of collinear mass 
singularities |J. The use of this theorem is normally based on two approximations. First, 
one picks out the leading-twist contribution and expresses the dimensionless cross section 
F(x,Q 2 ) ~ Q 2 a(x,Q 2 ) for a hard hadronic process as a convolution of pro cess- dependent 
coefficient functions C a and universal parton densities f a (Fig. 1), as follows 

F(x, Q 2 ) = C a (x, a 5 (/i 2 ), Q 2 /V) ® fa(x, v 2 ) + 0(A 2 QCD /Q 2 ) . (1) 

Here higher twists, i.e. contributions from multi-parton initial states (Fig. lb), are ne- 
glected, since they are suppressed by inverse powers of the hard scale Q 2 , and the leading- 
twist parton densities f a (x, /i 2 ) fulfil the renormalization group evolution equations 

- "j^lT " = Pab(a s (fi 2 ),x) ® / b (x,/i 2 ) , (2) 

where P a b{cxsi %) are the generalized splitting functions, whose Mellin transforms define the 
anomalous dimensions 7 a 6,jv 

7ab,N(a s ) = [ dx x N P ab (a s ,x) = P ab , N+ i(a s ) , (3) 
Jo 

N being the moment conjugate to x. 

Second, one considers the perturbative expansion for the kernels P ab of the yU 2 -evolution 

P ab (a s ,x) = £ (^^^(i) , (4) 

as well as the analogous expansion for the coefficients C a , and evaluates them to fixed order 
in as- Typically, QCD perturbation theory is at present under control up to two loops, 
that is, the terms P^ [|3] and P^ [^||| in the splitting functions (f|) are fully known, and 
the coefficient functions for most of the relevant processes at colliders are also known to 
the corresponding accuracy M. 

In the regime of low x, however, both the leading-twist approximation and the fixed- 
order truncation of evolution kernels and coefficient functions become critical. On the 
one hand, the leading-twist treatment violates the unitarity bound on the deep inelastic 
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cross section at asymptotic energies. This is related to the fact that, as the cross sec- 
tion increases, higher-twist contributions from parton recombination and rescattering [[7|,[| 
become corrections of relative order 0(1) in the formula ([!]). Therefore they cannot be 
neglected, and in fact they are essential for asymptotic unitarity to be restored. On the 
other hand, large high-energy contributions of the type a^\n m 1/x affect the perturbative 
expansions (d) (as well as the analogous expansions for the hard coefficients) to all orders in 
as- They may spoil the convergence of the perturbative series at small x, as higher powers 
of as associated with multiple hard-jet emission may be compensated by large enhancing 
factors in ln(l/x). These corrections have to be identified and summed to all perturbative 
orders. 

Throughout this paper we will not discuss the issue of unit arizat ion, and stick to the 
leading-twist framework. Phenomenological support for this attitude is provided by the 
absence of any signal of unitarity corrections in HERA data, as well as by estimates 
showing that the onset of such corrections is expected to be well beyond the range of present 
colliders. We will instead focus on the high-energy logarithmic effects in perturbation 
theory. Such effects arise from multiple radiation of space-like gluons in the t-channel of the 
process. Each two-gluon intermediate state (such as that in Fig. 2) does indeed contribute 
a logarithmic factor lnx for x — > 0, or equivalently, in the moment space of Eq. (^), a pole 
1/N for N —> 0. The resulting higher-order structure of the splitting functions P a b is 

Pab~ l \x) ~ \ [ln n_1 x + C(ln n - 2 x)] , x -»• . (5) 

Analogous terms show up in the coefficient functions C a , but here we will concentrate on 
the case of the splitting functions. The key approach to dealing with these potentially large 
corrections is based on perturbative resummation. One can consider improved perturba- 
tive expansions, which systematically sum classes of leading, next-to-leading, etc., small-x 
logarithms (or iV-poles) to all orders in as- For instance, in the moment space one has the 
following expansion for the anomalous dimensions 



1ab,N{as) = 



k=l L 



a s x ' 



(6) 



The coefficients (leading), (next-to- leading) , etc., define the logarithmic hierarchy 
at high energy (i.e., small x). Once these coefficients are known, they can be combined 
with expansions of the type @ (after subtracting the resummed logarithmic terms in 
order to avoid double counting), to obtain a prediction throughout the region of x where 
as ln(l/x)^l (or as/N&l), which is much larger than the domain asln(l/x) -C 1 where 
the as-expansions (§) are applicable. 

The outline of this paper is as follows. In Sect. 2 we review the present leading-order 
knowledge of QCD at small x, and discuss the role of sub-leading effects. In Sect. 3 we 



present the results of a study of next-to-leading corrections [|T0| , pT|| . Sect. 4 contains the 
summary and some prospects. 



2. Power counting at small x 

On the basis of the definition @, let us briefly go through a simple power counting at 
high energy. It is worth recalling first that no super-leading terms ag/N k , 2n > k > n, are 
present in the anomalous dimensions. This effect is associated with the coherence of eluon 
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radiation in space-like processes at small x [12-14|. Super- leading terms do contribute to 



exclusive final-state distributions, but they cancel out in inclusive quantities. The issue 
of the exclusive structure at small x has been studied in Ref. |l3j (some phenomenological 
investigations have been carried out in Ref. ]T5| ) , and will not be touched upon here. We 
will then be concerned with single-logarithmic corrections at small x. 

Note next that flavour non-singlet observables do not couple to pure-gluon intermediate 
states. Therefore they are always less singular than 1/x at small x, i.e., in the moment 
space, they are regular for N —>■ order-by-order in as- All the high-energy contributions 
dg/N k (n > k > 1) are thus associated with the flavour singlet sector. 



The singlet evolution equations are 

d ( f s \ _( 2N fl s m 2N flqg 



dln/i 2 V f g J V Igq Igg 




(7) 



where we have used standard definitions for parton densities and anomalous dimensions, 
and neglected regular terms at iV — » in the quark-quark entry of the matrix. 

It is a fundamental result which can be traced back to the work of Lipatov and collab- 
orators JT2] that only the gluon channel (lower entries to the anomalous dimension matrix 
in Eq. (|7])) contributes to the leading order (A^-type coefficients in Eq.@), and that the 
leading-order resummation can indeed be carried through by means of an integral equa- 
tion for the off-shell gluon Green function (BFKL equation). The resulting anomalous 
dimensions read as follows (as = Gaols I ^) 

J 99 ,nM = lQ) + 0(a s (^r) k ) , lgq!N (a s ) = g j 9g , N (*s) + 0(a s Q) k ) , (8) 
where the BFKL anomalous dimension j(as/N) is determined by the implicit equation 

1 = ^ XilQ)) , X(T) = 2^(1) - m - V(l - 7) , (9) 

ijj being the Euler ^-function. By solving Eq. (|) as a power series in the coupling constant 
one finds the leading (as/N) k contributions to the gluon anomalous dimension to all orders 
in as- The first perturbative terms are 

^)^ + ^(l) 4+ ^(f^ o ((f) , ) ■ m> 

£ being the Riemann ^-function (C(3) — 1.202, ^(5) ~ 1.037). Note that due to strong 
cancellations between real emission diagrams and virtual corrections the coefficients of a|, 
a s and a| in Eq. (|ToD vanish, and therefore the BFKL anomalous dimension departs from 



its one-loop approximation rather slowly. It is also worth recalling |16[ that by virtue of the 
collinear regularity of the BFKL integral kernel the leading-order anomalous dimensions 
are factorization scheme invariant, i.e., they do not depend on the explicit procedure 
to regularize and factorize collinear mass singularities. The physical features of the BFKL 
resummation may be schematically described considering the asymptotic behaviour of the 
gluon density f g . In the fixed-coupling limit we find for x — ► 

fg(x, Q 2 IQl) ~ exp (2y/a s In 1/x In Q 2 /Q%) (one-loop) (11) 
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f g (x,Q 2 /Ql) ~ x~ N (Q 2 /Ql) , iV~41n2a s (BFKL-resummed) (12) 

The all-order resummation of the perturbative iV-poles (|10|) builds up a branch point 
singularity at a value N of the moment proportional to as- As a result, the gluon density 
increases at small x like a power rather than like exp(^/ln(l/a;)). Also, the associated 

scaling violations are stronger in the resummed case ((Q 2 / Ql) 1 ^ 2 vs. exp[\Jhi(Q 2 / Qq)]) . 
The BFKL analysis predicts a growth of the anomalous dimension with respect to the 
fixed-order case, and eventually its saturation at the asymptotic value 1/2. 

The evaluation of small-x contributions to the anomalous dimension matrix (|7]) beyond 
the leading logarithmic approximation has been the object of many efforts over the past 
fifteen years. As in any perturbative expansion, one needs to know sub-leading corrections 
in order to be able to assess the stability of the expansion and set its limits of validity. In 
particular, one may wonder whether the asymptotic singularity (|i~2"|) may change beyond the 
leading order. Moreover, since running coupling effects mix with sub-leading corrections, 
it is essential to know the latter in order to be able to carry out the renormalization group 
analysis consistently at low x. 

It is worth noting that beyond the leading order quarks start to contribute on the same 
footing as gluons: corrections 0{as{ct s /N) k ) (B^-type coefficients in Eq.(|])) affect the 
lower as well as the upper entries to the anomalous dimension matrix in Eq. (ffh. This 
remark is relatively trivial, but it has been long overlooked in the literature, and most 
studies have focused on the pure-gauge sector. Let us give a perturbative example at fixed 
order. Consider the off-shell quark Green function (Fig. 3), and its cts-expansion. In the 
lowest non-trivial order 0(a s ), the remark above simply amounts to the well-known fact 
that the diagram in Fig. (3b) gives rise to a leading- order contribution (as/N) to the gluon 
anomalous dimension (times a term 0(as) in the coefficient function) when evaluated in 
the region of ordered momenta q 2 ^> k 2 ^> p 2 , but it contributes a next-to-leading- order 
correction ct 2 s /N to the quark anomalous dimension from the "disordered" configuration 
q 2 ~ k 2 ^> p 2 . A mechanism of this kind holds in higher-order diagrams as well. 

Actually, from a phenomeno logical viewpoint, one may argue that the knowledge of 
the next-to-leading quark anomalous dimensions is even more relevant than that of the 
corresponding corrections to the gluon anomalous dimensions. The reason for this is that 
the most accurate information on small-x QCD is coming from HERA data on deep inelastic 
structure functions, which couple to quarks directly, and to gluons via a C(as)-suppressed 
coefficient function. 

The next-to-leading contributions 0{as{cts/N) k ) to the quark anomalous dimensions 
(upper entries to the matrix in Eq. (|7])) have recently been computed to all orders in as 
fl0|]llll . In the next section we present the results of this calculation. The next-to-leading 
corrections to the gluon anomalous dimensions are in contrast still unknown at present. A 
calculational program is however being pursued by Fadin and Lipatov ]l7fl . 

3. Next-to-leading resummed results 

Unlike the leading-order analysis, the resummation of next-to-leading logarithms at 
small x is sensitive to the specific procedure one uses to regularize and factorize the collinear 
mass singularities arising from the low-momenta region. One then needs to develop a 
formalism in which high-energy resummation fl8|] is consistently matched with the all-order 
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collinear analysis. This can be accomplished [] l.| by virtue of a property of factorization at 
high energy, which holds for off-shell Green functions and supplements the one due to the 
renormalization group. 

We have performed the resummation of the next-to-leading corrections to the quark 
anomalous dimensions in two of the most commonly used factorization schemes, the DIS- 
and the MS-scheme. The DIS-scheme is defined by setting the coefficient functions for the 
structure function F 2 equal to unity in the quark channel and zero in the gluon channel, 
modulo an additional condition which is needed to uniquely fix the gluon density to all 
loops 0,11]. In this scheme the resummed result for the next-to- leading quark anomalous 
dimensions is parametrized as follows 



N 



l qg ,N{as) 



(13) 



where the anomalous dimension T^jvf can be expressed in a particularly compact form: 

-R(1n) ■ (14) 



' 19v 

, as m 2 + 3 lN -3j 2 N r 3 (l- 77 v)r 3 (l + 77 v) 



l qg ,N M = ^T R 



3 - 2 lN T(2 + 2 lN ) T(2 - 2 7 . 



N 



Here 7 tv denotes the BFKL anomalous dimension 
given by 



and the normalization factor is 



R{1n) 



T(l - In) x{1n) 
,T(1 + 7iV ) [-1nX'(1n)\ 



exp < 7jv0(1) + dl 



^)'{l) -^\l-~l) 



x{i) 



(15) 



X, x' being the characteristic function in Eq. (|S|) and its first derivative, respectively. The 
resummation of the contributions as{as/N) k to all orders in as is incorporated in Eq. (|14]) 
through the explicit ^-dependence and the as/iV-dependence of 7 at (known from the 
BFKL equation (0)). For instance, using the expansion (0), one can compute the first 
perturbative terms of the anomalous dimension (TIT 



ADIS) 



+ 




13 as 
If AT 



+ 



71 



-C(2) 



i?Z8-n C(2) + ?l C (3)-6C(4) 
81 18 SV ; 9 SV ' SV ; 




2.30 



Nj 



as\' 
NJ 

as_V 
NJ 



+ 



27 6 SV ; 3 SV . 



+ . . . 



+ 



,27 ,^ 



I 3 + 14.92 (^ * 



+ • 



asY 
NJ 



(16) 



Here the coefficients of the first two terms in the curly bracket agree with the known one- 
and two- loop anomalous dimensions in the DIS scheme P,|2"0| 
terms represent new sub-leading information at small x. 



whereas the higher-order 



Observe that Eq. (13) does not introduce any singularity in the N- moment space above 
the leading one (see Eq. (0))- This means that the position of the leading singularity is 
not changed by next-to-leading corrections in the quark sector. However, the approach to 
it is made faster, as a consequence of the positive sign of the corrections, and the resulting 
scaling violations are stronger. 
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It is worth noting that, as long as the gluon anomalous dimensions are unknown in next- 
to-leading order, the DIS-scheme may be convenient for phenomenological investigations 
of the structure function F 2 at small x, because by definition it decouples F 2 from gluons. 
The knowledge of the quark anomalous dimensions in this scheme (Eqs. (O)-(Il)) may 
therefore be particularly useful. 

As for the MS-scheme, the resummed result is expressed by relations of the type (p~3|) 
in terms of the analogous anomalous dimension 7q^jv (^sO- This is determined implicitly 



as a function of as and N by an algebraic equation (see Ref.|TT[ for full details), and its 
perturbative expansion reads 



+ 




- 

2tt 3 



! 535 14/05 
3 N 9 V N 



+ 



1.67^ + 1,6 



- + 2C(3) 
81 SV ; 



atsY 
N) 



+ 3.42 



N 



+ 5.51 



122 
243 



<xs\ 
N) 



(17) 



+ 



Note that the coefficients in Eq. fllTf ) are systematically larger than those in Eq. (|Ti 
This is related to the fact that the difference between the two anomalous dimensions is 
proportional to the MS-coefficients of the structure function F 2 , which are quite sizeable 
at small x [|Il[ . 



One may ask what the consequences of resummation are in terms of the stability of the 
logarithmic expansion at low x. One must first notice that the numerical implementation 
of QCD evolution based on resummed anomalous dimensions is not available yet, and more 
work is needed. Moreover, a fully consistent analysis to next-to-leading logarithmic order 
obviously requires also the computation of the gluon anomalous dimensions, which are still 
unknown to this accuracy. However, one can already make some remarks based on fixed- 
order features. First, the cancellations in orders 0(a s ), 0(a s ), 0(a s ) that we observe 



in the leading gluon anomalous dimensions (|T0| ) do no longer occur in the next-to-leading 
quark anomalous dimensions (|IBD (or (|T?D). Thus at large but finite energies (such as those 
at HERA) we may expect the latter to be of comparable importance to the former from 
the numerical point of view, in spite of their being formally sub-leading []1U|] . Second, a 
detailed analysis [fjl|] of the third-loop term of the anomalous dimension (|T7|) has recently 
been performed, which confirms this expectation. The authors of Ref. have considered 
the full two-loop evolution of parton densities (as in standard sets of distributions, like the 
MRS' p2|), and the evolution including the 0(a s ) term in eq. (|i~7]) and the 0(a s ) term in 
eq. (|K]). They have compared the predictions for F 2 in the two cases. They have done this 
both in the case of a flat input gluon distribution at low momenta (as in the set MRSDO) 
and a steep one (as in the set MRSD_). The outcome is that the impact of the higher-order 
corrections heavily depends on the choice of the input. If the input gluon distribution is 
flat, the effect may be sizeable already in the range of x- values accessible at HERA. This 
suggests a possible instability of the perturbative series, and calls for a careful evaluation of 
sub-leading contributions to all orders. On the other hand, a very small effect is observed 
in the case of a steep input gluon: this simply means that if the input is more singular than 
any rise that may ever be generated in QCD, then the input dominates the evolution, and 
the output at large momentum scales essentially reproduces what has been assumed as an 
input at low scales. 
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4. Conclusions 

QCD anomalous dimensions at small x are resummed to leading logarithmic accuracy 
by the BFKL equation. Sub-leading corrections to this approximation are needed for both 
theoretical and phenomenological reasons. In this paper we have reported on a study aimed 
at investigating such corrections. 

The analysis of the small-x region beyond the leading order brings in the novel feature 
of the interplay between high-energy logarithms and collinear mass singularities. To deal 
with this, a method has been set up in which the small-x resummation is carried out in a 
manner which can unambiguously be matched with the leading-twist collinear factorization 
to all orders. 

Sub-leading corrections to the BFKL analysis involve the quark and the pure-gauge 
sector on an equal footing. The main achievement of the work presented here is the resum- 
mation of the next-to-leading corrections to quark anomalous dimensions. Explicit results 
have been obtained in the MS and DIS factorization schemes. On the other hand, the 
corresponding contributions in the gluon sector are still unknown. 

Preliminary numerical studies indicate that the size of the next-to-leading effects dis- 
cussed in this paper may be large enough to affect the deep inelastic measurements at 
HERA. This suggests that further and more detailed numerical investigations should be 
pursued (including the full implementation of resummed results), as well as further theo- 
retical analyses on sub- leading contributions (in this respect, the calculation of the next- 
to-leading corrections to gluon anomalous dimensions represents a major challenge). Also, 
the procedure of matching mentioned above between the high-energy formalism and the 
collinear factorization should be used systematically to combine small-x resummed results 
with finite-x non-logarithmic contributions computed in fixed-order perturbation theory. 

Acknowledgments. The work presented in this paper has been carried out in collab- 
oration with S. Catani. Many useful discussions with M. Ciafaloni, G. Marchesini and B. 
Webber are gratefully acknowledged. 
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